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Standard methods for including electromagnetic interactions in lattice quantum chromodynamics 
calculations result in power-law finite-volume corrections to physical quantities. Removing these by 
extrapolation requires costly computations at multiple volumes. We introduce a photon mass to 
alternatively regulate the infrared, and rely on effective field theory to remove its unphysical effects. 
Electromagnetic modifications to the hadron spectrum are reliably estimated with a precision and 
cost comparable to conventional approaches that utilize multiple larger volumes. A significant 
overall cost advantage emerges when accounting for ensemble generation. The proposed method 
may benefit lattice calculations involving multiple charged hadrons, as well as quantum many-body 
computations with long-range Coulomb interactions. 
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Introduction — Approximately 95% of the visible 
mass of the Universe arises from the binding of quarks 
into nucleons by the strong interactions of quantum chro¬ 
modynamics (QCD). The relative mass difference be¬ 
tween the proton and neutron is approximately 0.07%, 
and is attributed to two sources of isospin symmetry 
breaking in the standard model, namely, differences in 
the down and up quark masses and their electric charges. 
Although these breaking effects are minute, they play an 
essential role in our understanding of the Universe. For 
example, the primordial abundance of light nuclear ele¬ 
ments in the early Universe is exquisitely sensitive to the 
excess mass of the neutron compared to the proton mm- 

Lattice QCD (LQCD) provides a first-principles 
approach for determining isospin-breaking effects in 
hadronic and nuclear processes. There are a handful of 
LQCD calculations of the strong contribution to the nu¬ 
cleon mass splitting w and a comparable number that 
determine the electromagnetic corrections lillSHTO]. One 
impressive calculation includes both sources of isospin 
breaking simultaneously and yields, among other quanti¬ 
ties, a postdiction for the nucleon isospin splitting with 
^ 5(7 statistical significance [S] . There exists an alternate 
means for determining the electromagnetic self-energy 
of the nucleon, from the Cottingham formula [nH2D], 
which makes use of experimental cross sections as in¬ 
put to dispersion integrals. However, the uncertainty at¬ 
tained with this method ^11 - 05] is not yet competitive 
with the LQCD calculations. 

Although inclusion of electromagnetism in LQCD is 
theoretically straightforward [211 US] , it presents practi¬ 
cal challenges due to the long-range nature of the elec¬ 
tromagnetic (QED) interactions. Specifically, such inter¬ 
actions give rise to power-law finite-volume (FV) correc¬ 


tions, and their removal via extrapolation requires com¬ 
putationally demanding simulations performed at multi¬ 
ple volumes. An analytic understanding of the power-law 
FV effects within such setups |8l [26U^ has enabled re¬ 
liable FV extrapolations of the single hadron spectrum. 

Despite the successful application of present tech¬ 
niques, there are a number of reasons for considering new 
methods. Control over FV modifications to light nuclear 
binding energies seems to require particularly large vol¬ 
umes |26] . There are quantities in addition to the spec¬ 
trum for which a precise knowledge of the QED modi¬ 
fications is needed, for example, corrections to hadronic 
matrix elements [22] and charged particle scattering [3D] , 
both of which suffer from infrared (IR) challenges. LQCD 
calculations are performed with multiple ultraviolet (UV) 
regulators, providing valuable cross-checks on the con¬ 
tinuum extrapolation of many important quantities Ell- 
Multiple IR regulators can do the same for LQCD calcu¬ 
lations that include QED, but to date, only a few other 
formulations have been considered [331 - I33| . Of those, 
only one is constructed with a local quantum field the¬ 
ory [33] [31]. Finally, computationally efficient means of 
accounting for IR effects are always desirable, not just for 
lattice QCD-I-QED, but anywhere long-range Coulomb 
interactions are present (see, e.g., Ref. [33]). 

Motivated by these considerations, we demonstrate 
the viability of an alternative IR regulator for lattice 
QCD-I-QED simulations: namely, the introduction of a 
photon mass m-y. Although a photon mass term mani¬ 
festly violates gauge invariance, it maintains locality and 
its effects on hadronic quantities can be reliably quan¬ 
tified and accounted for within an effective field theory 
(EFT) framework. The introduction of a new scale, m-y, 
implies an additional extrapolation within our approach. 


2 




0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 


FIG. 1. Zero-mode adjusted (filled) and unadjusted (open) ef¬ 
fective mass difference for the kaon splitting (m-i/fhr: = 1/14 
and Lja = 24). The diagonal grid lines have slope 2x\ the 
red and blue points correspond to different sinks. The gray 
bands correspond to uncertainties on the extracted value for 

AMeff. 


With the aid of analytic formulas, however, we demon¬ 
strate that for the spectrum, such extrapolations can be 
performed at a single volume and yield results that are 
consistent with conventional approaches. In the remain¬ 
ing sections, we present the salient features of our calcu¬ 
lation. 

Analytic considerations — In continuum Euclidean 
spacetime, the gauge fixed action for the massive pho¬ 
ton is given by 

■^7 = 

where — dyA^-, throughout this study, we 

work in Landau gauge, corresponding to the limit ^ 0. 

An Abelian theory, such as QED, with a massive vector 
gauge field is still perturbatively renormalizable. This 
well known result follows from the fact that it is possible 
to find a Becchi-Rouet-Stora-Tyutin (BRST) transforma¬ 
tion that leaves the Lagrangian invariant up to a total 
divergence |86) . The BRST symmetry is preserved if one 
uses a gauge invariant UV cutoff m, such as a space- 
time lattice; thus, the renormalizability follows from the 
power-counting theorems for a lattice regularization |38j . 

We consider three forms of corrections to correlators 
and hadron mass differences at leading order in the fine- 
structure constant a = e^/(47r). These corrections arise 
from either the zero mode contribution to the partition 
function, the presence of a finite photon mass, or FV 
effects. The analytic forms of these corrections are de¬ 
termined from an EFT for hadrons of mass M {M = mn, 
nip, nix+, and m^o) and charge Q; the naive expansion 
is in m-y/M (i.e. Ayv = M) [SS]. The EFT is a gen¬ 
eralization of nonrelativistic QED m for hadrons that 
includes a photon mass term, and additional operators 
that are not constrained by gauge invariance. 

(1) Zero mode: For sufficiently small m^, the zero mode 
of the temporal photon field appearing in Eq. must be 


FIG. 2. QEDj,^ induced mass differences, extrapolated to 
infinite volume (taking Kl = 2). 


treated nonperturbatively |41] . In this regime, the two- 
point function for single hadrons has the form 

C'(t) = , (2) 


where Z is an overlap factor; the zero-mode contribu¬ 
tion appears as a; = (47ra(5^)/(2m^L^T), and vanishes 
as T —> oo at fixed L and 

(2) Photon mass: The hadron’s electromagnetic mass 
shift can be determined as a function of photon mass, or¬ 
der by order in an expansion in powers of m^/M. With 
the electromagnetic mass written as M{a,my), we de¬ 
fine the mass shift AyM{a,mj) = AI{a,my) — M{a,0), 
which is UV finite. These IR shifts are given by 

^ _ « Q.) A ,3) 

The leading-order (LO) expression is nonanalytic in 
the squared photon mass, whereas the next-to-leading- 
order (NLO) expression is analytic but arises from both 
loops and local contributions [15]; the next-to-next-to- 
leading-order (N^LO) correction is of order A^M^ = 
0{m^/M‘^). The latter two orders are accompanied by 
coefficients not fixed by the hadron charge. 

(3) Finite volume: The effects of FV can similarly be 
calculated using a nonrelativistive QED approach. This 
is a finite photon mass generalization of that pursued by 
Refs. [26l|27]. The FV corrections to the electromagnetic 
mass are written as 5LM{a,m^,L) = M{a,m^,L) — 
M{a,m-y,oo), and for charged hadrons are given up to 
NLO by 


5lM^^ = 2'KaQ^my 




[niyLY 


TTl 

S^j^nlo ^ [2Ii/2(to^L) +l3/2{ni^L)] ,(4) 


M 


where 


d^n{z) 


1 

2"+W5r(n) 




(5) 
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TABLE I. QEDj,^ induced mass splittings, extrapolated to 

1/ —>■ 00 (m-y = 0, — 0). 


splitting 

Kl 

xVdof 

AM/M X 10^ 

p — n 

1 

0.07/2 

0.73(05) 


2 

0.03/1 

0.70(13) 

4- 

1 

o 

1 

0.29/2 

3.71(06) 


2 

0.17/1 

3.68(20) 


TABLE 11. QEDJ^^ induced 
rrij —>■ 0 {L/a = 24, Kl = 1) 

mass 

splittings, extrapolated to 

splitting 

m-f/m-s: range 


X^/dof 

AM/M X 10^ 

p — n 

1/14 - 1 

2 

0.09/5 

0.79(06) 


1/4 - 1/2 

1 

0.06/2 

0.81(08) 

K+ -K^ 

1/14 - 1 

2 

0.42/5 

3.77(06) 


1/4 - 1/2 

1 

0.12/2 

3.79(06) 


and v G By contrast, the leading nonvanishing cor¬ 
rection for neutral baryons [mesons] appears at N^LO 
[next-to-next-to-next-leading-order {N^LO)]. Because 
the zero mode of the temporal photon is treated exactly 
in Eq.[^ the FV corrections are calculated with this mode 
removed-a manifestation of which is the subtracted term 
appearing at LO. 

Lattice parameters and ensembles Electro- 
quenched numerical calculations of the hadron spectrum 
were performed using a modified version of the Chroma 
software suite [43]. Studies were performed using dy¬ 
namical SU(3) flavor symmetric isotropic QCD gauge 
field configurations generated using a tadpole-improved 
Liischer-Weisz gauge action and clover fermion action. 
The configurations correspond to a single lattice spac¬ 
ing a = 0.1453(16) fm, three spatial extents: L ^ 3.48 
fm, 4.64 fm and 6.96 fm, and temporal extents T > L. 
The pion (kaon) and nucleon masses in physical units are 
m^=fhK = 807.0(9.1) MeV and m„ = 1.634(18) GeV, 
respectively. This choice of masses ensures that the only 
appreciable FV corrections to hadron masses are those 
arising from QED effects. The QCD ensembles used 
in this work comprise 956 {L/a = 24,T/a = 48), 515 
(L/a = 32, T/a = 48) and 342 {L/a = 48, T/a = 64) 
configurations and are a subset of those described in 
Ref. [44]; further details regarding the ensembles, lattice 
action and parameters can be found there. 

The uncorrelated photon field configurations were 
generated using two different lattice actions: a con¬ 
ventional massless Coulomb gauge-fixed action with the 
zero mode removed [TTl [M] (QEDj.^) |3S|, and a 
naive lattice discretized form of Eq. (QED^), where 
the derivatives are replaced by finite differences. Note 
that in Euclidean space, Landau gauge is a complete 
gauge-fixing condition, and therefore in the latter case, 
the path integration over nonzero modes is well de¬ 
fined in the -G 0 limit. The photon mass val¬ 


ues considered in this work are given by m^/rfiT^ € 
[1/14,1/7,1/4,1/3, 5/12,1/2, 7/12,1]. In both cases, re¬ 
sults were obtained by computing correlation functions 
on QCD-I-QED gauge conhgurations generated by post- 
multiplying each QCD configuration by a single , 

where Qu = 2/3, Qd = Qs = —1/3. In the elec- 
troquenched approximation with SU(3) flavor symme¬ 
try, isospin splittings have missing contributions that are 
O(a^), and therefore negligible for this study. 

In the electroquenched theory, the hne structure cou¬ 
pling does not renormalize and therefore we take it to 
be equal to its experimental value a~^ = 137.036..., 
measured in the Thomson limit. The presence of elec¬ 
tromagnetic interactions demands renormalization of the 
valence bare quark masses m^, however. Since our lat¬ 
tice regulator breaks chiral symmetry, this leads to an 
additive shift in the quark mass. We tune the valence 
quark masses so that, in the presence of electromagnetic 
interactions, the neutral qq meson mass niqq obtained 
from the connected part of the qq correlation function is 
sufficiently close to the pion (kaon) mass For our 

electroquenched calculation, this choice of renormaliza¬ 
tion is robust but the quark mass renormalization in the 
full QCD-fQED does not allow for a unique separation of 
the QED and QCD effects [15]. All measurements were 
performed using valence quark masses am„ = —0.25501 
and amd = arris = —0.24750 (the QCD bare quark mass 
is arriq = —0.2450); the resulting mistuning for the charge 
neutral mesons was Amqq/ffiTr ^ 0.1% for all values of 
m-y/ffiTc < 1, where Arriqq = rriqq — ffiT^. 

The mistuning from strong isospin breaking can be es¬ 
timated using chiral symmetry. For the kaon, one finds 


friK 


1 Arriuu - Arridd 

2 iflK 


< 0.0004, 


( 6 ) 


while the nucleon correction is given by 


Arrin-p 2{Amdd - Ato„„) ml 

z — (^d—u z 


47r/^m„ 


(7) 


We can estimate the parameter ad-u from the LQCD 
determination of the m^-contribution to the nucleon 
mass splitting HHH] and find Amp_„/m„ 0.0002. In 
both cases, mistuning is a potentially sizable correction 
to our results, which affects both the QED-j.^ and QEDj^^ 
determinations. Although a precise quark mass tuning is 
required for practical applications, it is not needed in the 
present proof-of-principle study [T7] . 

Analysis and results — Shell-shell and shell-point 
correlation functions were estimated using a single mea¬ 
surement per configuration, with a randomly chosen 
spacetime source location. Following Ref. |9], we aver¬ 
age observables over +e and —e on each configuration in 
order to exactly cancel off the 0(e) contributions to sta¬ 
tistical noise. Mass differences due to electromagnetic ef¬ 
fects can be determined from the late-time dependence of 
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FIG. 3. Volume adjusted (Kl = 1) QCDm induced mass differences, extrapolated to m-y —>■ 0. Fits were performed using 
L/a = 24 data at all (right) and only the middle four (left) values of rrij. 


the single hadron correlation functions C^{t) and C^{t), 
by studying the plateau region of an effective mass dif¬ 
ference = M^{t) — By exploiting 

the correlations between A and B, we are able extract a 
clear signal for the mass difference. For the nucleons, we 
consider a generalized effective mass formula of the form: 


.^eff.exp (f ) 


1 

a 


C{t a) 


2xt + xa , 


( 8 ) 


which neglects the backward propagation of states on a 
lattice of finite temporal extent T. For mesons, we ac¬ 
count for the backward propagating state by considering 
a generalized effective mass formula of the form: 




- cosh ^ 
a 


- ^h(T,a) _|_ ^h{T,-a)- 

_ 


xT,{9) 


where h{T,a) = xa{a — T 2r) -h log[C'(r -h a)/C{T)]. 
Both formulas treat the zero mode of the temporal pho¬ 
ton field appearing in Eq. nonperturbatively (for neu¬ 
tral hadrons cc = 0 and these expressions reduce to their 
conventional forms). Although this contribution is negli¬ 
gible compared to the hadron masses, for the lattice pa¬ 
rameters considered it can be comparable in magnitude 
to the mass differences we wish to extract. Fig. pro¬ 
vides an explicit example of the behavior of AMeff(r) for 
the kaon mass splitting, computed both with and without 
the zero-mode contribution accounted for. 

Mass differences were determined for all volumes and 
photon masses via a correlated constant least-squares fit 
to AMeff in the plateau region, as demonstrated in Fig.j^ 
An analogous determination from exponential fits to a 
ratio of correlation functions yielded consistent results. 
Systematic uncertainties were estimated by varying the 
region over which fits were performed, and all uncertain¬ 
ties were added in quadrature. The extracted mass shifts 
were subsequently extrapolated to vanishing photon mass 


and/or the infinite-volume limit using the fit formula: 

AM {a, L, m~y) = AM {a) + ^ AjM^ 

k=0 

Kl 

{aym^yL), (10) 

where Kj and Kl indicate the order of each extrapola¬ 
tion. In the case of mass splittings, an appropriate linear 
combination of mass shift formulas was used. Note that 
for the extrapolations, = 0; the appropri¬ 

ate EV formulas for 5lM^'‘^^ retain T dependence, and 
may be found in Ref. [5]. 

We carry out two independent analyses to test the vi¬ 
ability of our proposal: (1) an infinite-volume extrapo¬ 
lation of QEDy^ induced mass differences, as is conven¬ 
tionally performed, and (2) an —)■ 0 extrapolation of 

QEDj\^ induced mass differences using data at a single 
EV, but after having first removed the lowest order EV 
contributions 5lM. Both types of extrapolation were 
performed using Eq. |10[ noting that many of the lowest- 
order contributions are fixed by theory. Results for the 
first analysis, using all three volumes, are provided in Ta¬ 
ble]^ and representative fits are shown in Eig. Results 
for the second analysis on the smallest volume are pro¬ 
vided in Table |ll] for comparison, and shown in Eig. 
Analogous QED^ extrapolations, performed at each of 
the three volumes, are summarized in Fig.j^and are con¬ 
sistent not only with each other, but also the QED^^ 
extrapolations. In all cases, we find that the numerical 
and theoretical mass corrections are in excellent agree¬ 
ment down to at least m-yL 1. 

The most computationally demanding part of our cal¬ 
culation involves multiple inversions of the Dirac oper¬ 
ator. Assuming, conservatively, a linear scaling with 
spacetime volume, the total inversion cost for L/a = 32 
is 515/956 X (32/24)^ 1.3 times greater than that of 

L/a = 24. By comparison, the L/a = 48 inversion cost 
is 4^ 3.8 times greater. The L/a = 24 extrapolations us¬ 
ing my/friTy e [1/4,1/2] data, provided in Table [Hj are 
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FIG. 4. QEDjyj^ extrapolations (points) performed indepen¬ 
dently at each volume. Systematic errors are estimated by 
considering fits over multiple ranges of and orders < 3 
(summarized by the histograms). QEDj.^ extrapolations re¬ 
sults are summarized by gray horizontal bands. 

consistent with those using all values of m^. The results 
are also consistent with the extrapolation using 

three volumes, provided in Table [Ij but required only 4/5 
the computational cost. We therefore conclude that for 
the same precision and accuracy, the numerical cost of 
our QED^ calculation of the mass splittings is compa¬ 
rable to or less than that of QEDj^j;^. 

Conclusion — This work demonstrates that it is pos¬ 
sible to reliably estimate infinite-volume hadron mass 
differences induced by electromagnetism on a single lat¬ 
tice volume with QED^. Conservatively, the pionless 
EFT employed in this work is valid for <C 2m,r and 
Wti-E 4. Provided these inequalities are satisfied, the 
analytic expressions obtained for the mass shift are valid 
up to 0{m^/M^, a'^), and are independent of the pion 
mass; from our numerics, it appears that this order is 
sufficient to obtain reliable extrapolations of the mass 
shifts in the regime m-y/rriT^ ^ 1 and m-yL > 1. 

On preexisting lattice configurations, and for equal 
computational cost, we obtain an equally precise uncer¬ 
tainty in extrapolated differences as compared to the tra¬ 
ditional method. This cost comparison does not account 
for the significant overhead of generating the configura¬ 
tions in the first place. The results of our analysis pave 
the way for a more complete treatment of QED correc¬ 
tions using this approach. When considering more in¬ 
volved LQCD calculations, such as charged-particle scat¬ 
tering [30] , our method provides a mass gap to produce a 
photon, thus increasing the range of energy for which the 
standard Liischer method [ISl SS] for obtaining the scat¬ 
tering phase shift can be employed. It will be interesting 
to explore these types of calculations, and also to use our 
method with chiral fermions, which do not suffer from 
additive quark mass renormalization. Finally, it would 
be interesting to see if our method of screened interac¬ 
tions coupled with analytic extrapolation techniques is of 
benefit to quantum many-body calculations. 
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